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Introduction 

Over the last few years, there has been a revival of interest in VTOL aircraft capable 
of operating in fixed wing as well as rotary wing mode. High speed rotorcraft designs, 
such as the tilting rotor configuration, pose an entirely new problem in the rotary wing 
field. The design goals for this class of aircraft include low downwash velocity in hover, 
good low speed maneuverability and cruise speeds of 350 - 500 knots 1 . Several new 
concepts 2 ' 5 have recently been proposed to meet these design goals. Extensive research 
performed in this field have led to the XV- 15 research aircraft and ultimately to the 
production of the V-22 Osprey tilting rotor for the US Navy. 

The combined requirements of efficient high speed performance of a fixed wing 
aircraft and good helicopter-like hover characteristics complicates the design process of 
tilting high speed proprotor aircraft. It is necessary to maintain good aerodynamic 
efficiency in high speed axial flight without degrading hover efficiency. This often leads 
to conflicting design requirements. For example, improved efficiency in high speed 
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cruise demands high drag divergence Mach numbers which are normally associated with 
thin airfoils. This however, reduces the hover figure of merit by reducing Cj/c. 

Therefore, to maintain the required thrust ceiling in hover, the rotor solidity has to 
increase. Also as the forward speed increases, helical tip Mach number limitations 
caused by whirl flutter, require a reduction in the rotor rotational velocity. Introducing 
blade sweep can alleviate this problem by reducing the effective chordwise Mach 
number, which allows for higher speeds, without reducing the rotor RPM. Therefore the 
proper design of proprotor blades capable of achieving the design objectives must 
consider the right combination of airfoil thickness and blade sweep in addition to other 
aerodynamic variables such as planform and twist. 

Several studies have been performed 6 - 9 to study design trade offs between the two 
flight modes. For example, Johnson et al. 6 performed a detailed study on the 
performance, maneuverability and stability of high speed tilting proprotor aircraft, 
including the XV- 15 and V-22. Liu and McVeigh 7 recently studied the use of highly 
swept rotor blades for high speed tilt rotor use. However, formal optimization techniques 
were not used in these studies. Recently an effort was initiated by Chattopadhyay and 
Narayan 8 - 9 to develop formal multidisciplinary optimization procedures for the design of 
civil high speed tilting proprotor blades. The propulsive efficiency in axial flight was 
maximized with constraints on the figure of merit in hover, aeroelastic stability in cruise 
and other aerodynamic and structural design criteria. The purpose of the present paper is 
to formulate the optimum design problem of high speed proprotors using multiobjective 
optimization techniques with the integration of the necessary disciplines. 

Problem Definition 

f / y 

A integrated, multiobjective optimization procedure is developed for the design of 
j high speed proprotors with the coupling of aerodynamic, dynamic, aeroelastic and 
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structural criteria. The objectives are to maximize propulsive efficiency in high speed 
cruise and rotor figure of merit in hover. Constraints are imposed on rotor blade 
aeroelastic stability in cruise and on total blade weight. Two different multiobjective 
formulation procedures, the Min EP and the K-S function approaches are used to 
formulate the two-objective optimization problem. 

Aerodynamic Model 

The rotor studied is a wind tunnel model an existing advanced technology proprotor, 
which is a three bladed rotor with a rigid hub. Cubic variations are assumed for the chord 
(c) and twist (9) distributions to model the blade aerodynamics, 

c(y) o c 0 + ci(y - 0.75) + c 2 (y - 0.75) 2 + c 3 (y - 0.75) 3 (1) 

0(y) = 9 0 + 0i (y - 0.75) + 0 2 (y - 0.75) 2 + 0 3 (y - 0.75) 3 (2) 

where y denotes the nondimensional blade radius. Note that c D represents the chord and 
0 O the twist at the 75 percent blade radius, respectively. A quadratic lifting line is used 
and is defined as follows. 

x = f(y) = eiy + £ 2 y 2 (3) 

where Ei, e 2 are constants that determine the curvature, and are defined such that 

I £i I £ Ci (4) 

where £i (i = 1, 2) are prescribed bound for the curvature parameters. These bounds 
allow for either forward or backward in-plane curvature. When ej and e 2 are equal to 
zero the lifting line will be a straight line. The blade sweep, based upon this lifting line 
distribution, assumes the following form 
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where A(y) is the sweep distribution, in degrees, defined to be positive aft of the straight 
lifting line. 


Structural Model 

The structural model used for the problem consists of a simple two-celled box beam 
as shown in Fig. 1. The beam is considered to be the principal load carrying member of 
the proprotor and the stiffness contributions from the honeycomb and the nonstructural or 
tuning weights are placed at the blade tip and are distributed along the planform. The 
total blade weight comprises the weight of the box beam, the skin, the honeycomb and 
the nonstructural weights. The wall thicknesses of the box beam are assumed to vary in 
proportion to the chord distribution. 



Figure 1 Double-celled box beam configuration 
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Optimization 

Objective — Eunctions: The multiobjective optimization procedure is used to 

simultaneously maximize the rotor propulsive efficiency in high speed cruise and the 
hover figure of merit. 

Design Variable ! : Both aerodynamic and structural design variables are used. The 
aerodynamic design variables include chord, twist and sweep distribution coefficients 
(Eqns. 1, 2 and 5). The structural design variables comprise roots values of the wall 
thicknesses of the two-cell box beam and the magnitudes of the nonstructural weights, at 
the tip, distributed spanwise. 

. Constraints : To avoid the possible occurrences of air and/or ground resonance associated 
with a soft inplane rotor, it is important to maintain the value of the lowest natural 
frequency in hover, fj, above 1/rev. Therefore the following constraint is imposed. 

(i) fi > 1/iev (g) 

It is important to impose aeroelastic stability constraints to prevent any degradation of the 
rotor stability in high speed cruise. This is all the more important when the blade mass 

and stiffness are altered during optimization. The stability constraints are expressed as 
follows. 

(ii) a k < -D k k = 1, 2, ... , K (7) 

Where K represents the total number of modes considered, and a k is the real pan of the 
stability root. The quantity \) k denotes the minimum allowable blade damping and is 
defined to be a small positive number. To avoid incorporation of weight penalties, after 
optimization, the total blade weight is constrained as follows. 

(iv) W < Wy 


( 8 ) 
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Multiobjective Optimization 


A typical optimization problem involving multiple objective functions can be 
mathematically posed as follows. 


Minimize Ffcftjin) 


Subject to 

gj(<t>n) ^ 0 

constraints) 


k = 1, 2, ... , NOBJ (objective functions) 
n= 1,2, ... ,NDV 

j = 1,2, ... , NCON (inequality 


4>nL ^ ^n ^ (^nu (side constraints) 

where NOBJ denotes the number of objective functions, NDV is the number of design 
variables and NCON is the total number of constraints. The subscripts L and U denote 
lower and upper bounds, respectively, on the design variable (|> n . A description of the 
multicriteria design objective formulation follows. 

This study examines three multiobjective function formulation techniques that are 
less judgmental than the Pareto based weighting factors and are therefore more suited to 
large scale, highly nonlinear optimization problems that are associated with rotary wing 
design. The two multiobjective function techniques used are the Minimum Sum Beta 
(Min ZP) and the Kreisselmeier-Steinhauser (K-S) function approaches. A description of 
these methods follows. 

Minim um .Sum Beta (Min Zpi : This method was first used by Weller at al. 10 to 
formulate a two objective function rotor vibration problem. Using these technique, 
pseudo-design variables that represent tolerances of the individual objective functions 
from prescribed tolerances are introduced. The objective function, F](<J)), is then defined 
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as a linear combination of these tolerances of each objective function to their specified 
target values as follows 

NOBJ 

Fi(<D) = ^p k ( 9) 

k=l 

where Pk are pseudo design variables with properties such that the original objective 
functions Fk remain within a Pk tolerance of some prescribed values. This requirement 
introduces new constraints of the following form. 

£ p k k= 1,2,..., NOBJ (10) 

F k 

The quantities Fk are the prescribed target values of the individual objective functions Fk. 
Using the above formulation, as the objective function, and correspondingly the values of 
Pk, are reduced to zero, the values of the individual objective functions Fk are driven to 
their prescribed values, Fk- The design variables for the Min Lp formulation comprise 
the original set of design variables and the pseudo design variables, Pk- 

Kreisselmeier-Stei nhauser (K-S) function : This technique was first utilized by Sobieski 
et al. 11 at the NASA Langley Research Center. The first step in formulating the objective 
function in this approach involves transformation of the original objective functions into 
reduced objective functions These reduced objective functions take the form 

* Fk(O) 

Fk(0) = ' L °' gmax " 0 k=1 > > N OBJ (11) 

where Fk 0 represents the value of Fk calculated at the beginning of each iteration. The 

quantity g max is the value of the largest constraint corresponding to the design variable 
vector O and is held to be constant for each iteration. These reduced objective functions 
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are analogous to the previous constraints, and therefore a new constraint vector g2 m (<I > ) 
(m = 1 , 2 , ... , M) is introduced, where again M = NCON + NOBJ. The new objective 
function to be minimized is then defined, using the K-S function as follows: 

M 

F 2 (<x>) = fmax +-ln ^ e p(gm (<),) " fma ^ (12) 

m=l 

where f max is the largest constraint corresponding to the new constraint vector, g2 m (0), 
and in general is not equal to gmax- The multiplier p can be considered analogous to a 
draw-down factor where p controls the distance from the surface of the K-S objective 
function to the surface of the maximum constraint function. When p is large the K-S 
function will closely follow the surface of the largest constraint function. When p is 
small, the K-S function will include contributions from all violated constraints. The 
design variable vector <J> is identical to that used in the Min Ef} approach. 

Analysis 

Dynamic, Aerodynamic and Aeroelastic Analyses The aerodynamic, dynamic and 
aeroelastic analysis of the high speed proprotor is performed using the code 
CAMRAD/JA^. The code has the capability of analyzing both helicopter and tilting 
rotor aircraft. Wind tunnel trim options are used as the reference blade is a wind tunnel 
blade model. In cruise, the blade is trimmed to specific rotor lift and drag coefficients 
using the rotor collective and cyclic pitch angles. A prescribed wake model, as 
implemented in CAMRAD/JA, is used to model the aerodynamics in hover and the rotor 
is trimmed to a specific value of the coefficient of power. In axial flight, the components 
of the induced velocity are negligible compared to the high forward speed of the rotor. 
Therefore, uniform inflow conditions are used to model the aerodynamics in this case. 
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The aeroelastic stability analysis in cruise is analyzed with a constant coefficient 
approximation. 

Structural Analysis: The detailed structural analysis of the rotor blade is performed 
based upon the two-celled trapezoidal box beam using an inhouse code that was recently 
developed specifically for this application. 


Optimization Implementation 

The optimization is performed by using the program CONMIN 14 . The program uses 
the method of feasible directions to solve nonlinear constrained optimization problems. 
To reduce the computational effort, an approximate analysis technique is used to compute 
the objective function and constraint values during iterations within the optimizer. For 
this problem the two-point exponential hybrid approximation technique ^ is used. This 
technique takes its name from the fact that the exponent used in the expansion is based 
upon gradient information from the previous design point. The technique is formulated 
as follows. 


£(<*>) = F(<D 0 ) 



3F(O n ) 

&t>n 


( 13 ) 
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8F((Di) 


Pn = 


log 


s \ 

d<J > 0 


J) 


+ 1.0 


n J 


( 14 ) 


The quantity <J>j refers to the design variable vector from the previous iteration and the 
quantity Oq denotes the current design vector. The exponent p n can be considered as a 
“goodness of fit” parameter, which explicitly determines the trade-offs between 
traditional and reciprocal Taylor series based expansions (also known as a hybrid 
approximation technique). Details of this method can be found in Ref. 15. 

Results and Discussions 


A wind tunnel model of an existing high speed proprotor is used as a baseline design. 
The optimization for this problem is performed with a cruise velocity of 400 knots and a 
rotational velocity of £l = 375 RPM (tip speed of 491 ft/s) in axial flight. The operating 
condition is 20,000 feet above sea level. In hover, a rotational velocity of Cl = 570 RPM 
(tip speed of 746 ft/s) is used at sea-level conditions. The high forward flight speed of 
400 knots represents the target cruise value for high-speed rotorcraft. The tip speed is 
reduced in the airplane mode so that the helical tip Mach number stays below Mjd (the 
drag divergence Mach number). The rotor RPM in cruise (375) is selected after 
performing a parametric study on the effect of forward speed and rotor RPM on 
propulsive efficiency. A value of Cj/g = 0.08 is used to trim the blade in forward flight, 
and a value of Cp/a = 0.0131 is used to trim the blade in hover. The blade radius is 12.5 
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feet, and the blade is discretized into 10 segments (NSHG = 10). For the Min 10 

approach 24 design variables are used (including the pseudo-design variables), and in the 
K-S function approach 22 design variables are used. 

Some results, obtained to date, are summarized in Table 1 and Figs 2 and 3. Table 1 
presents a summary of preliminary optimization results. From Table 1 and Fig. 2 it can 
be seen that it is possible to obtain substantial increases in both the hover figure of merit, 
(21.7 - 28.8 percent), and the axial propulsive efficiency, ti ax (24.6 - 41.3 percent) using 
both multiobjective formulation techniques. It is of interest to note that the mean chord 
(and correspondingly the blade solidity) is increased by 71 percent and 40 percent in the 
K-S function and Min 10 approaches, respectively from the baseline value. Two possible 
explanations exist for this large increases in the rotor solidity. First, in order to satisfy the 
frequency constraint, the root chord is significantly increased to make the stiffnesses 
larger, which in turn increases the solidity (see Fig. 3). Secondly, since the hover figure 
of merit is being maximized, o is being increased to increase the thrust margin of the 
rotor in hover. 

Based on the previous experience, the above problem will be formulated with 
constraints on rotor solidity. The final paper will present results of the integrated 
aerodynamic/dynamic/aeroelastic optimization problem of high speed proprotors with 
additional design constraints. Design trade-off studies will also be performed by varying 

the flight conditions and the results of corresponding optimum blade designs will be 
presented. 
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Table 1 Summary of Preliminary Optimization Results 



Reference 

blade 

Bounds 

lower upper 

Optimum 

Min £p K-S 

Objective 

Functions 

FOM 

0.662 

- 

- 

0.853 

0.936 

'Hax 

0.647 

- 

- 

0.787 

0.805 

Constraints 

W (lb) 

194 

- 

194 

167 

173 

fi (per rev) 

0.812 

1.00 

- 

1.01 

1.34 

ai 

0.096 

- 

-0.001 

-0.040 

-1.529 

a 2 

0.096 

- 

-0.001 

-0.040 

-1.529 


-0.697 

- 

-0.001 

-0.732 

-0.169 

04 

-0.697 

- 

-0.001 

-0.732 

-0.169 

05 

-2.431 

- 

-0.001 

-2.443 

-2.502 

06 

-0.170 

- 

-0.001 

-0.265 

-0.073 

Pi 

0.150 

0.001 

0.200 

0.006 

- 

P2 

0.150 

0.001 

0.200 

0.010 

- 

Mean chord 

Ce (ft) 

1.48 

- 

- 

2.07 

2.52 

Solidity 

a 

0.113 

- 

- 

0.158 

0.193 

Trim 

C j/a 

0.110 

- 

- 

0.117 

0.116 


□ Reference 

□ Min LP 
■ K-S 
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